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Flame Propagation is used as a prototypical example of ex- 
panding fronts that wrinkle without limit in radial geometries 
but reach a simple shape in channel geometry. We show that 
the relevant scaling laws that govern the radial growth can 
be inferred once the simpler channel geometry is understood 
in detail. In radial geometries (in contrast to channel geome- 
tries) the effect of external noise is crucial in accelerating and 
wrinkling the fronts. Nevertheless, once the interrelations be- 
tween system size, velocity of propagation and noise level are 
understood in channel geometry, the scaling laws for radial 
growth follow. 
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manner. We will exemplify the details of such a trans- 
lation in the context of premixed flames that exist as 
self sustaining fronts of exothermic chemical reactions in 
gaseous combustion. But our contention is that similar 
ideas should be fruitful also in other contexts of unstable 
front propagation. Needless to say, there are aspects of 
the front dynamics and statistics in the radial geometry 
that cannot be explained from observations of fronts in a 
channel geometry; examples of such aspects are discussed 
at the end of this Letter. 

Mathematically our example is described [8] by an 
equation of motion for the angle- dependent modulus of 
the radius vector of the flame front, R(0, t): 
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The main idea of this Letter is that in order to de- 
rive scaling laws for unstable front propagation in radial 
geometry, it is useful to study noisy propagation in chan- 
nel geometries, in which the noiseless dynamics results 
usually in simple shapes of the advancing fronts [1] . Fa- 
mous examples of such situations are Laplacian growth 
patterns and unstable flame fronts. Laplacian growth in 
a channel geometry results in a single finger whose width 
is uniquely determined by the channel width L, the ve- 
locity of the tip v and the surface tension [2]. Lapla- 
cian growth in radial geometries results in convoluted 
and complex structures whose full characterization still 
eludes repeated theoretical attacks [2-5]. Flame prop- 
agation offers a similar situation; in channel geometries 
the flame front looks like a giant cusp whose velocity v is 
constant [6] . In radial geometries the flame front acceler- 
ates with time, and while propagating it wrinkles by the 
addition of hierarchies of cusps of all sizes [7-12]. 

The understanding of radial geometries requires con- 
trol of the effect of noise on the unstable dynamics of 
propagation. It is particularly difficult to achieve such a 
control in radial geometries due to the vagueness of the 
distinction between external noise and noisy initial condi- 
tions. Channel geometries are simpler when they exhibit 
a stable solution for growth in the noiseless limit. One 
can then study the effects of external noise in such geome- 
tries without any ambiguity. If one finds rules to trans- 
late the resulting understanding of the effects of noise 
in channel growth to radial geometries, one can derive 
the scaling laws in the later situation in a satisfactory 
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Here 0< 8 < 2tt is an angle and the constants Ub,DM 
and 7 are the front velocity for an ideal cylindrical front, 
the Markstein diffusivity and the thermal expansion co- 
efficient respectively. Ro(t) is the mean radius of the 
propagating flame: 
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The functional I(R) is best represented in terms of its 
Fourier decomposition. Its Fourier component is 
where Rk is the Fourier component of R. Simulations of 
this equation, as well as experiments in the parameter 
regime for which this equation is purportedly relevant, 
indicate that for large times Rq grows as a power in time 



Ra(t) — (const + t) 



13 



(3) 



with (3 > 1, and that the width of the interface W grows 
with Rq as 



w(t) ~ R (ty 
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with \ < 1. Note that this exponent may differ signifi- 
cantly from the exponent which characterizes the scaling 
of correlation functions [13]. This phenomenon is com- 
mon to cases in which large features dominate the width 
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of the graph. In a previous publication [11] we showed 
that (i and \ satisfy the scaling relation (3 = 1/x- In this 
Letter we determine the numerical value of these expo- 
nents from the analysis of noisy channel growth. 

In channels there is a natural lengthscale, the width 
L of the channel. The translation of channel results to 
radial geometry will be based on the identification in the 
latter context of the time dependent scale C(t) that plays 
the role of L in the former. To do this we need first to 
review briefly the main pertinent results for noisy channel 
growth. In channel geometry the equation of motion is 
written in terms of the position h(x, t) of the flame front 
above the x-axis. After appropriate rescalings [6] it reads: 



v~L», ^« 0.35 ±0.03 



(9) 



dh(x,t) _ 1 
dt ~2 



dh{x,t) 



dx 



d 2 h{x,t) 

dx 2 



+ I{h(x,t)} + 1. 



(5) 



It is convenient to rescale the domain size further to < 
9 < 27r, and to change variables to u(9,t) = dh(9,t)/d9. 
In terms of this function wc find 
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dt 
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where L = L/2n. In noiseless conditions this equation 
admits exact solutions that are represented in terms of N 
poles whose position Zj (t) = Xj (t) + iyj (t) in the complex 
plane is time dependent: 
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u(9,t) = t^y^cot 



9- Zj (t) 



c.c. 
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The steady state for channel propagation is unique and 
linearly stable; it consists of N(L) poles which are aligned 
on one line parallel to the imaginary axis. The geometric 
appearance of the flame front is a giant cusp, analogous 
to the single finger in the case of Laplacian growth in a 
channel. The height of the cusp is proportional to L, and 
the propagation velocity is a constant of the motion. The 
number of poles in the giant cusp is linear in L, 
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The introduction of additive random noise to the dy- 
namics changes the picture qualitatively. It is conve- 
nient to add noise to the equation of motion in Fourier 
representation by adding a white noise rjk for every k 
mode. The noise correlation function satisfies the rela- 
tion < r) k {t)r) k >(t') >= 5(k + k')5(t - t')2f/L. The noise 
in our simulations is taken from a flat distribution in the 
interval [— y/2f /L, y/2f / L] ; this guarantees that when 
the system size changes, the typical noise per unit length 
of the flame front remains constant. It was shown [12] 
that for moderate but fixed noise levels the average ve- 
locity v of the front increases with L as a power law. In 
our present simulations we found 



For a fixed system size L the velocity has also a power 
law dependence on the level of the noise, but with a 
much smaller exponent: v ~ , £ « 0.02. These 
results were understood theoretically by analysing the 
noisy creation of new poles that interact with the poles 
defining the giant cusp [12]. It was shown that the sys- 
tem is always linearly stable or marginally stable for the 
introduction of new poles, but with increasing L all the 
eigenvalues (which are negative for a positive measure of 
L values) all decrease in absolute magnitude like l/L 2 , 
making the system more and more susceptible to noisy 
perturbations [12]. The picture used remains valid as 
long as the poles that are introduced by the noisy per- 
turbation do not destroy the identity of the giant cusp. 
Indeed, the numerical simulations show that in the pres- 
ence of moderate noise the additional poles appear as 
smaller cusps that are constantly running towards the 
giant cusp. Our point here, is not to predict the nu- 
merical values of the scaling exponents in the channel 
(this was done in ref. [12], but to use them to predict the 
scaling exponents characterizing the acceleration and the 
geometry of the flame front in radial geometry. 

Superficially it seems that in radial geometry the 
growth pattern is qualitatively different. In fact, close 
observation of the growth patterns (see Fig.l) shows that 
most of the time there exist some big cusps that attract 
other smaller cusps, but that every now and then "new" 
big cusps form and begin to act as local absorbers of small 
cusps that appear randomly. The understanding of this 
phenomenon gives the clue how to translate results from 
channels to radial growth. 




FIG. 1. Simulations of the outward propagating flame 
front. Note that there is a wide distribution of cusp sizes. 
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In the radial geometry we non-dimensionalize Equation 
(2) by rescaling R — > r, Rq — > r and t-»r, and finally 
rewriting the equation in terms of u = dr/d9: 

To complete this equation we need a second one for r (t) . 
This equation is obtained by averaging (2) over the angles 
(i.e. operating on the equation with ^ f Q * dd): 

f\He + i. (ii) 

dr 2r 2 2irJ Q K 1 

These equations again admit exact solutions in terms of 
poles, of the form of Eq.(7). It is easy write down the 
equations of motion of the poles and check that the poles 
are attractive along the real direction (which means phys- 
ically that they are attracted along the angular coordi- 
nate) but they are repulsive along the imaginary direc- 
tion, which is associated with the radial coordinate. If 
it were not for the stretching that is caused by the in- 
crease of the radius (and with it the perimeter), all the 
poles would have coalesced into one giant cusp. Thus we 
have a competition between pole attraction and stretch- 
ing. Since the attraction decreases with the distance be- 
tween the poles in the angular directions, there is always 
an initial critical length scale above which poles cannot 
coalesce their real coordinates when time progresses. 

Suppose now that noise adds new poles to the system. 
The poles do not necessarily merge their real positions 
with existing cusps. If we have a large cusp made from 
the merging of the real coordinates x c of N c poles, we 
want to know whether a nearby pole with real coordi- 
nate x\ will merge with this large cusp. The answer will 
depend of course on the distance D = r \x c — xi\. A 
direct calculation [11], using the equation of motions for 
the poles shows that there exists a critical length C(ro) 
such that if D > C(r ) the single pole never merges with 
the giant cusp. The result of the calculation is that 

C~ry p . (12) 

Note that a failure of a single pole to be attracted to a 
large cusp means that tip-splitting has occurred. This is 
the exact analog of tip-splitting in Laplacian growth. 

It is now time to relate the channel and radial geome- 
tries. We identify the typical scale in the radial geometry 
as C <~ W ~ . On the one hand this leads to the scaling 
relation (3 = On the other hand we use the result 

established in a channel, (9), with this identification of a 
scale, and find f = . Comparing with (3) we find: 

'-ra- <13) 

This result leads us to expect two dynamical regimes 
for our problem. Starting from smooth initial conditions, 



in relatively short times the roughness exponent remains 
close to unity. This is mainly since the typical scale C is 
not relevant yet, and most of the poles that are generated 
by noise merge into a few larger cusps. In later times the 
roughening exponent settles at its asymptotic value, and 
all the asymptotic scaling relations used above become 
valid. We thus expect (3 to decrease from 1/(1 — /i) to an 
asymptotic value determined byx=l//?in(13): 

13 = 1 + u « 1.35 ±0.03 . (14) 

The expected value of \ 1S thus x = 0.74 ± 0.03. 
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FIG. 2. Panel a: a logarithmic plot of the velocity ver- 
sus time for a radially evolving system. The parameters of 
the simulation are: / = 1CP 8 , 7 = 0.8, v = 1. Panel 
b:Logarithmic plot of the width of the flame front as a func- 
tion of the mean radius. 

We tested these predictions in numerical simulations. 
We integrated Eq.(6), and in Fig. 2 we display the re- 
sults for the growth velocity as a function of time. After a 
limited domain of exponential growth we observe a con- 
tinuous reduction of the time dependent exponent. In 
the initial region we get (3 = 1.65 ±0.1 while in the final 
decade of the temporal range we find (3 = 1.35 ± 0.1. Wc 
consider this a good agreement with (14). A second im- 
portant test is provided by measuring the width of the 
system as a function of the radius, see Fig. 2b. Again we 
observe a cross-over related to the initial dynamics; In the 
last temporal decade the exponent settles at \ = 0.75±.l. 
We conclude that at times large enough to observe the 
asymptotics our predictions are verified. 

Finally, we stress some differences between radial and 
channel geometries. Fronts in a channel exhibit mainly 
one giant cusp which is only marginally disturbed by the 
small cusps that are introduced by noise. In the radial 
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geometry, as can be concluded from the discussion above, 
there exist at any time cusps of all sizes from the smallest 
to the largest. This broad distribution of cusps (and 
scales) must influence correlation function in ways that 
differ qualitatively from correlation functions computed 
in channel geometries. 
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mate F'(y) — ^2e =y P(£)- The second derivative yields 
F"{y) « —P(y). Thus the structure function is deter- 
mined by the scale distribution of cusps, and if the latter 
is a power law, this should be seen in the second deriva- 
tive of F(y) as demonstrated in Fig. 3. The conclusion 
of this analysis is that the radial case exhibits a scal- 
ing function that characterizes the distributions of cusps, 
P(£) « £- a . 

In summary, we demonstrated that it is possible to 
use information about noisy channel dynamics to pre- 
dict nontrivial features of the radial evolution, such as 
the acceleration and roughening exponents. It would be 
worthwhile to examine similar ideas in the context of 
Laplacian growth patterns. 
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FIG. 3. Panel a: a logarithmic plot of the correlation func- 
tion of the interface. Panel b: Second derivative of the corre- 
lation function of the interface. 

To make the point clear we exhibit in Fig. 3a the struc- 
ture function 



F(y) = y/(\R(x + y)-R( X )\*) 
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computed for a typical radial front, with x = R9. To 
stress the scaling region we exhibit the second deriva- 
tive of this function in Fig. 3b. The low end of the graph 
can be fitted well by a power law y~ a with a w 0.6. 
This indicates that F(y) w Ay + By 2 ~ a . In a chan- 
nel geometry we get entirely different structure functions 
that do not exhibit such scaling functions at all. The 
way to understand this behavior in the radial geometry 
is to consider a distribution of cusps that remain distinct 
from each other but whose scales are distributed accord- 
ing to some distribution P{£). For each of these cusps 
there is a contribution to the correlation function of the 
form f(y,£) ~ £g(y/£) where g(x) is a scaling function, 
g(x) w x for x < 1 and g(x) « constant for x > 1. The 
total correlation can be estimated (when the poles are 
distinct) as 
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F(y)^P(e)£g{y/e) . 
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The first derivative leaves us with ^ t P{£)g' '(y / '£), and 
using the fact that g' vanishes for x > 1 we esti- 



4 



